CO oxidation on a single Pd atom supported on magnesia 



O 
O 
(N 



(N 



C/3 



X3 

o 



> 
o 

o 



I 

o 
o 



X 



S. Abbet,^ U. Heiz,^ H. Hakkinen'^ and U. Landman'^ 
^Universite de Lausanne, Institut de Physique de la Matiere Condensee, CH-1015 Lausanne, Switzerland 
^University of Ulm, Institute of Surface Chemistry and Catalysis, 89069 Ulm, Germany 
^School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430 
(to be published in PRL 18 June 2001) 

The oxidation of CO on single Pd atoms anchored to MgO(lOO) surface oxygen vacancies is 
studied with temperature-programmed-reaction mass-spectrometry and infrared spectroscopy. In 
one-heating-cycle experiments CO2, formed from O2 and CO preadsorbed at 90 K, is detected at 260 
K and 500 K. Ab-initio simulations suggest two reaction routes, with Pd(CO)202 and Pd(C03)C0 
found as precursors for the low and high temperature channels, respectively. Both reactions result 
in annealing of the vacancy and induce migration and coalescence of the remaining Pd-CO to form 
larger clusters. 

PACS: 82.65.-hr, 68.43.Bc, 68.47.Jn 



Molecular-scale understanding of the energetics and 
mechanisms of catalytic reactions could open new av- 
enues to the design of catalysts with specific functions 
To this aim model catalysts are used |^ to ex- 
tract information on reactivities at conditions relevant 
to catalysis [||-D. A most fruitful approach to gain in- 
sights into the reaction mechanisms of catalytic processes 
is the combination of experiments and high-level atomic 
scale simulations. Indeed, such recent joint studies re- 
vealed the reaction mechanism of the oxidation of CO on 
Aug clusters and the polymerization of acetylene on 
Pd„ (n < 30) clusters deposited on MgO(lOO) films §. 

Here we report on studies of the oxidation of CO 
by a model catalyst consisting of single Pd atoms an- 
chored on oxygen surface vacancies (F-centers, PCs) 
of a MgO(lOO) film. After low-temperature (90 K) 
coadsorption of O2 and CO the formation of CO2 was 
detected by temperature-programmed-reaction (TPR) 
mass-spectrometry in single-heating-cycle experiments. 
Ab-initio density-functional simulations were performed 
to identify relevant molecular precursors as well as to 
study the CO2 formation mechanisms. The determi- 
nation of two initial molecular complexes adsorbed on 
MgO(lOO), Pd(CO)202 and Pd(C03)C0, is supported 
by good agreement between the measured and calculated 
CO vibrational frequencies. The former complex is in- 
volved in the formation of CO2 at 260 K, and the de- 
composition of the carbonate complex leads to CO2 des- 
orption at 500 K. Both reaction routes induce annealing 
of the surface PC and migration of the remaining Pd-CO 
unit to form larger Pd clusters. 

Model Catalyst Preparation and Experimental Results. 
Size-selected atomic Pd cations were deposited on an in- 
situ prepared MgO(lOO) thin film |]. The low-kinctic- 
energy (0.2 eV) deposition of only a 0.45 x 10~^ of a 
monolayer (ML) of Pd at a substrate temperature of 90 
K reduces greatly the migration of the Pd atoms bound 



to the PCs (whose concentration is 10 ~^ ML, i.e., roughly 
twice that of the adsorbed Pd atoms). Indeed, a recent 
comparison of ab-initio calculations and FTIR studies of 
CO adsorbed on the supported Pd atoms provided clear 
evidence for single Pd atoms bound to the PCs of the 
MgO(lOO) support. The observation of two different vi- 
brational bands at 90 K suggests the presence of at least 
two different CO molecules adsorbed on monodispersed 
Pd atoms (see ref. Q and below). A recent study on 
the cyclotrimerization of acetylene also revealed that the 
PC- trapped Pd atoms are stable up to 300 K 

First, we verified that the clean MgO(lOO) thin films 
are inert for the oxidation reaction; i.e., no CO2 was 
formed in a one-heating-cycle experiment after adsorb- 
ing O2 and CO or vice versa |10|. When Pd atoms are 
trapped on the PCs, preadsorption of oxygen and subse- 
quent saturation by CO leads to the formation of carbon 
dioxide, with desorption peaks at 260 K and at around 
500 K (Fig. 1). The existence of two desorption peaks 
suggests the presence of two different reaction mecha- 
nisms. Note that when CO is preadsorbed prior to O2 
the oxidation reaction is suppressed, indicating CO poi- 
soning. 

Information pertaining to the mechanism of the CO 
oxidation on PC-trapped Pd atoms was obtained by mea- 
suring the CO vibrational bands during reaction (See in- 
sets to Fig. 1). At 95 K three features in the FTIR 
spectrum of ^^CO are observed. The infrared absorption 
at 2125 cm~^ originates from CO adsorbed on extended 
defect sites on the MgO(lOO) thin films The broad 
band with a peak at around 2045 cm^^ and a shoulder 
at 2005 cm^^ indicates adsorption of at least two CO 
molecules. Heating to 165 K results in a band at 2005 
cm~^ with a shoulder at 2035 cm~^. The decreased inten- 
sity of the high-frequency band suggests partial desorp- 
tion of the high frequency CO molecule prior to the oxida- 
tion reaction. At 250 K, close to the temperature of maxi- 
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mum CO2 desorption, the shoulder at 2035 cm~^ almost 
disappears, indicating oxidation of the ^'^CO molecule 
with a vibrational frequency of 2035 cm^^. Further heat- 
ing to 410 K results in disappearance of the band at 2005 
cm~^. This correlates with the formation of CO2 and the 
observed desorption of molecular CO (not shown) . In ad- 
dition, a new vibrational band with a frequency of 1830 
cm^^ appears between 250 and 300 K. The disappear- 
ance of this band above 600 K correlates with complete 
CO oxidation and molecular CO desorption. 

Theoretical Methodology and Results. The calculations 
were performed using the Born-Oppenheimer (BO) local- 
spin-density (LSD) molecular dynamics method (BO- 
LSD-MD) |12| w ith the generalized gradient approxima- 
tion (GGA) |13| and employing norm-conserving non- 
local scalar-relativistic [Q (for the Pd atom) pseu- 
dopotentials [|l5|. Such calculations yield accurate re- 
sults pertaining to geometries, electronic structure, and 
charging effects of various neutral and charged coinage 
metal clusters |0 and nanostructures jl^]. The mag- 
nesia surface was modeled by a finite region ("clus- 
ter" ) of atoms, whose valence electrons are treated fully 
quantum- mechanically (using the BO-LSD-MD), embed- 
ded in a large (2000 charges) point-charge lattice, as de- 
scribed in our study of the Aus /MgO model catalyst [|j . 

A single Pd atom binds strongly to the oxygen va- 
cancy (binding energy of 3.31 eV), with a slight amount 
of charge (0.15 e) transferred to the adsorbed atom. In 
comparison, the binding energy of Pd atoms to terrace 
oxygen sites is only 1.16 eV. The enhanced binding to 
the FCs is also reflected in the corresponding bonding 
lengths of 1.65 A and 2.17 A for MgO(FC)-Pd and MgO- 
Pd, respectively |l8[ . 

Binding of two CO molecules saturates the MgO(FC)- 
Pd system; occupying the MgO(FC)-Pd system with 
three CO molecules leads to spontaneous (barrierless) 
desorption of one of the molecules. In the most stable 
configuration the two CO molecules are inequivalent; one 
CO binds on-top and the second adsorbs on the side of 
the Pd-atom [9b] (this top-side geometry is similar to 
that shown in Fig. 2a but without the O2), and the total 
binding energy of the two CO molecules is 1.62 eV. An 
alternative symmetric adsorption configuration, with the 
two CO molecules adsorbing on opposing sides of the Pd 
atom, is less stable by 0.61 eV than the top-side one [9b]. 

To study the oxidation mechanisms of CO on 
MgO(FC)-Pd the system was optimized first with coad- 
sorbed O2 and two CO molecules. Two stable geomet- 
ric arrangements were found, with the most stable one 
shown in Fig. 2a where the CO molecules bind in a 
top-side configuration and the O2 is adsorbed parallel 
to the surface on the other side of the Pd atom. This 
configuration (with spin 5 = 0) is 0.90 eV more stable 
than an alternative one {S — 1) where the O2 is bound 
on-top of the Pd atom and the two CO molecules oc- 
cupy the side positions (not shown). The preadsorbed 



O2 molecule enhances slightly the adsorption energy of 
the two CO molecules (1.78 eV) compared to the case 
without preadsorbed oxygen (1.62 eV). The adsorbed O2 
molecular bond is stretched and activated (1.46 A com- 
pared to the calculated gas-phase value of 1.25 A). In ad- 
dition, we found a stable carbonate complex Pd(C03)C0 
(Fig. 3a), whose binding energy is 4.08 eV larger than 
the aforementioned Pd(CO)202 complex. 

Reaction mechanisms. Two reaction mechanisms are 
proposed corresponding to the two CO2 peaks observed 
experimentally (Fig. 1). At low temperatures the two 
relevant precursors are shown in Fig. 2a and Fig. 3a. 
The existence of the Pd(CO)202 complex (Fig. 2a) is 
supported by the agreement between the calculated and 
measured vibrational frequencies of the two inequiva- 
lently adsorbed CO molecules (see Table I) . The absence 
of a clear FTIR signal corresponding to the carbonate 
species can originate from the adsorption geometry (see 
Fig. 3a) where the carbonate is bound at the side of the 
Pd atom and the CO3 plane is only slightly tilted away 
from the surface, resulting in vanishingiy small normal 
dynamic dipole components; however, the frequency of 
the side-bonded CO of this complex (2020 cm~^) lies in 
the experimentally observed vibrational band. 

Corresponding to the 260 K CO2 desorption peak we 
propose the following reaction mechanism. First, in a 
competitive process, CO desorbs or is oxidized upon heat- 
ing. The theoretically estimated activation energies of 
the two processes are 0.89 eV for desorption (Table I) and 
0.84 for oxidation, obtained from a series of constrained 
energy-minimizations, where the top-CO molecule ap- 
proaches the closest O atom of the O2 molecule. The 
transition state leading to CO2 formation is shown in 
Fig. 2b. After finding the transition state we performed 
a microcanonical MD simulation to study the reaction 
dynamics. The desorbing CO2 molecule (Fig. 2c) car- 
ries away the major part (~2 eV) of the reaction heat of 
about 2.2 eV, partitioned as 0.1 eV, 0.1 eV, and 1.8 eV 
into the translational, rotational and vibrational degrees 
of freedom, respectively, i.e., the desorbing molecule is vi- 
brationally " hot" . The remaining O atom of the complex 
fills the O vacancy under the adsorbed Pd atom (Fig. 2c), 
releasing 2.8 eV (not shown in Fig. 2d). Concomitantly 
with the annealing of the O vacancy the binding energy 
of the Pd atom is largely reduced from 3.31 to 1.16 eV, 
as discussed above. Note also that the binding energy 
of the remaining CO molecule to the Pd atom increases 
(the calculated binding energy of a CO molecule to a Pd 
atom adsorbed on the terrace of MgO(lOO) is 2.29 eV). 
Consequently, migration of the Pd-CO unit, leading to 
formation of larger clusters, becomes energetically feasi- 
ble (the calculated diffusion barrier for a Pd atom on the 
MgO(lOO) terrace is 0.43 eV). Indeed, the observed vi- 
brational band with a frequency of 1870 cm^^ appearing 
between 250 K and 300 K (Fig. 1) is in close agreement 
with the calculated frequency of a CO molecule bridge- 
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bonded to an adsorbed Pd2 dimer; for a Pd2 adsorbed 
on a MgO(lOO) surface the calculated frequency is 1836 
cm~^, and for a Pd2 adsorbed at an FC it is 1877 cm~^ 
(Table I). Desorption of the adsorbed CO molecule oc- 
curs at temperatures below 600 K. The calculated bind- 
ing energy of the bridge-bonded CO to the dimer is 1.95 
eV when it is adsorbed at an FC and 2.87 cV when ad- 
sorbed on a terrace, indicating that the Pd2 dimer (with a 
bridge-bonded CO) is most likely to be located at an FC, 
as otherwise an even higher CO desorption temperature 
is expected. Such FCs are available since the density of 
the oxygen vacancies in our experiment is roughly twice 
that of the Pd atoms. 

Formation of CO2 at higher temperatures (correspond- 
ing to desorption around 50GK, Fig. 1) involves decom- 
position of the Pd(C03)C0 carbonate complex (Fig. 3a). 
This mechanism is observed in molecular dynamics sim- 
ulations where the temperature is controlled to 500 K by 
Langevin dynamics. After the transition state (see Fig. 
3b), overcoming an energy barrier of about 1 eV (see 
Fig. 3d), CO2 leaves the complex parallel to the MgO 
surface with a total kinetic energy of about 0.25 eV, dis- 
tributed approximately as 0.1 eV and 0.15 eV between 
the translational and vibrational modes, respectively, and 
with a vanishing rotational component. The remaining 
O atom fills the O vacancy (Fig. 3c) as found also for 
the lower-temperature CO2 formation mechanism. The 
total cxothermicity of this process is 1.8 eV (sec the sharp 
drop in Fig. 3d for t > 210 fs). As in the low-temperature 
mechanism, the remaining Pd-CO can migrate and coa- 
lesce to larger clusters. 

Finally, we remark that the relative abundance of the 
two surface precursors, Pd(CO)202 and Pd(C03)C0, 
underlying the above reaction mechanisms may be domi- 
nated by kinetic factors; e.g., formation of the carbonate 
complex, although energetically favorable, could be hin- 
dered by the requirement that the CO molecule will ap- 
proach the preadsorbed side-02 rather than the Pd atom, 
as well as by a significant reaction barrier. Further stud- 
ies of the low temperature formation of the carbonate 
complex, as well as efforts to identify such species (e.g. 
via HREELS), are warranted. 
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TABLE I. Calculated CO binding energies (Eb) and ^^CO 
vibrational frequencies (t*-'*''), compared to the experimental 
frequencies {lo"^'' , scaled for ^^CO from the values shown in 
Fig. 1 for ^"^CO). The calculated dissociation energy, equilib- 
rium bond length, and harmonic frequency of the gas-phase 
^^CO molecule are 11.06 eV, 1.141 A, and 2140 cm"\ respec- 
tively, compared to the experimental values of 11.09 eV, 1.128 
A, and 2170 cm . The CO vibrational frequencies are deter- 
mined through molecular dynamics simulations stretching the 
equilibrium CO bond by 1% and observing the dynamics over 
a few harmonic vibration periods. All the calculated lj values 
include a correction factor of 2170/2140. The error estimate 
for the caJcnlated fre(iU(>iici(>s is cm V 
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Eb (eV) 


LO^'' {an-') 


LO^'^P (cm-^) 


MgO(FC)-Pd-(CO)2-02 
MgO(FC)-Pd-(C03)-CO 
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MgO-Pd2-CO 


1.78^ 

1.95 
2.87 


2019 / 2088" 
1687 / 2020° 
1877'* 

isse'* 


2050 / 2080 

1870 
1870 



''per 2 CO molecules 
''side- / top-CO 

'^maximum of the CO3 frequencies / side-CO 
"^bridge-bonded CO 



FIG. 1. A TPR spectrum showing formation of CO2 on 
a MgO(FC)-Pd sample after preadsorption of O2 and satura- 
tion with CO at 90 K. The insets show the infrared spectra 
of adsorbed *^C0 after heating the sample to the indicated 
temperatures. All the spectra were recorded at 90 K. 



FIG. 2. (color) (a) 

Optimized structure of the MgO(FC)-Pd-(CO)202 complex. 

(b,c) selected configurations, and (d) the potential energy vs 
time, recorded in an ab-initio MD simulation where CO2 is 
formed from the complex shown in (a). The simulation starts 
from the transition state shown in (b). The potential en- 
ergy of the transition state is 0.84 eV above the optimized 
configuration shown in (a), (c) A snapshot at 210 fs, where 
the formed CO2 is desorbing and the remaining O atom from 
O2 molecule is moving towards the F-center. The pertinent 
structural parameters in (a) are: top-CO: d(C-0) = 1.149 A, 
d(Pd-C)=2.10 A, /(Pd-C-0)= 178°; side-CO: d(C-0)=1.164 
A, d(Pd-C)=1.91 A, Z(Pd-C-0)=166°; O2: d(0-0)=1.46 A, 
d(O-Pd)=2.09 A and 2.31 A. The Pd atom is black, Mg and O 
ions of the substrate are blue and red, respectively, adsorbed 
O2 is yellow, carbons are grey and their respective oxygens 
are purple and green. 

FIG. 3. (color) (a) Optimized struc- 

ture of the MgO(FC)-Pd-(C03)CO complex. (b,c) selected 
configurations, and (d) the potential energy vs time recorded 
in an ab-initio MD simulation where the carbonate complex 
is decomposing at around 500 K. Snapshots from this sim- 
ulation are shown in (b) at 160 fs (just after the transition 
state) and (c) at 240 fs. The PC is filled by an O atom in (c) . 
The pertinent structural parameters in (a) are: the side CO: 
d(C-O)=1.160 A, d(Pd-C)=1.90 A, Z(Pd-C-0)=166°; CO3: 
d(Pd-O)=2.07 A, d(C-0)=1.45 A, 1.23A, 1.29 A with the 
long C-0 bond forming to the Pd-bound oxygen. Colors as 
in Fig. 2, with the carbonate oxygens shown in yellow. 
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